First observation of 7/(1405) decays into / (980)7r° 
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The decays J/ip — > 77r + 7r _ 7r° and J ftp — > 77r°7r 7r° are analyzed using a sample of 225 million 
J/tp events collected with the BESIII detector. The decay of ?7(1405) — s> /o(980)7r° with a large 
isospin violation is observed for the first time. The width of the /o(980) observed in the dipion mass 
spectra is anomalously narrower than the world average. Decay rates for three-pion decays of the 
77' are also measured precisely. 

PACS numbers: 14.40. Be, 12.38.Qk, 13.25.Gv 
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A state near 1440 MeV/c 2 was discovered in 
pp annihilation at rest decaying to ry7r + 7r~ and 
was subsequently observed decaying to KKir 
Considerable theoretical and experimental efforts have 
since been devoted to understand its nature. It was 
proposed to be a candidate for a pseudo-scalar glueball @, 
|3|; the measured mass, however, is much lower than that 
obtained from lattice QCD calculations, which is above 
2 GeV/c 2 jij. Later, experiments produced evidence that 
this state was really two different pseudo-scalar states, 
the 77(1405) and the 77(1475). The former has large 
couplings to ao(980)7r and KKir, while the latter mainly 
couples to K*K. A detailed review of the experimental 
situation can be found in Ref. 

The nature of the well-established light scalars /o(980) 
and do (980) is also a matter of controversy. It is not clear 
whether they belong to the light scalar meson nonet or 
are examples of mesons beyond the naive quark model 
(eg . tetra-quark states, hybrids or KK molecules) 0- 
Hl| . The possibility of mixing between the /o(980) and 
a°(980) was suggested long ago, and its measurement 
sheds light on the nature of these two resonances p^ - [l7| . 

The thrcc-pion decays of the 77' have garnered attention 
because their branching ratios (Br) can probe isospin 
breaking [TH, [l9|. The ratios of the branching ratios 



(r± = Br{rj — > tt + it ir a )/Br(?]' — > tt + tt 77) and ro = 
Br{rj -)■ 3ir )/Br(r]' -> ir°ir°ri)) are related to the 
strange quark mass and SU(3) breaking [l8j . 

In this letter, we present the results of a study of 
J/tp — > 77T + 7r~7r° and Jjip — > 77r 7r°7r°. A single 
structure around 1.4 GeV/c 2 in the 7r + 7r~7r° (7r 7r 7r°) 
mass spectrum is observed, associated with a narrow 
structure around 980 MeV/c 2 in the 7r + 7r - (7r 7r°) mass 
spectrum. This analysis is based on a sample of (225.2 ± 
2.8) x 10 6 J/ip events 1201 accumulated in the Beijing 
Spectrometer (BESIII) [2lJ operating at the Beijing 
Electron-Positron Collider (BEPCII) [22|. 

BEPCII is a double-ring e + e~ collider designed 
to provide e + e~ collisions with a peak luminosity 
of 10 33 cm _2 s _1 at a beam current of 0.93 A. 
The cylindrical core of the BESIII detector consists 
of a helium-based main drift chamber (MDC), a 
plastic scintillator time-of- flight system (TOF), and a 
CsI(Tl) electromagnetic calorimeter (EMC), which are 
all enclosed in a superconducting solenoidal magnet 
providing a 1.0 T magnetic field. The solenoid is 
supported by an octagonal flux-return yoke with resistive 
plate counter muon identifier modules interleaved with 
steel. The acceptance of charged particles and photons 
is 93% over 47r stereo angle, and the charged-particle 
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momentum and photon energy resolutions at 1 GeV are 
0.5% and 2.5%, respectively. The BESIII detector is 
modeled with a Monte Carlo (MC) simulation based on 
GEANT4 [23l.[2^. 

The charged-particle tracks in the polar angle range 
| cos 6>| < 0.93 are reconstructed from hits in the MDC. 
Tracks that extrapolate to be within 20 cm of the 
interaction point in the beam direction and 2 cm in 
the plane perpendicular to the beam are selected. The 
TOF and dE/dx information are combined to form 
particle identification confidence levels for the 7r, K, 
and p hypotheses; each track is assigned to the particle 
type that corresponds to the hypothesis with the highest 
confidence level. Photon candidates are required to 
have at least 25 MeV and 50 McV of energy in the 
EMC regions |cos0| < 0.8 and 0.86 < |cos0| < 0.92, 
respectively, and be separated from all charged tracks by 
more than 10°. 

For J/ijj — > 77r + 7T~7r°, the candidate events 
are required to have two oppositely charged tracks 
identified as pions and at least three photon candidates. 
A four-constraint(4C) energy-momentum conserving 
kinematic fit is performed to the 7777r + 7r~ hypothesis, 
and xlc < 30 is required. For events with more 
than three photon candidates, the combination with the 
smallest \ 2 is retained. To reject possible background 
events with two or four photons in the final state, the 
4C-fit probability for an assignment of the J/tp — > 



~7r 77 and the J/tp 



~tt 7777 channels. The 



7T° candidates are selected by requiring \M 1 



TO^O < 



0.015 GcV/c 2 . Events with \M 



< 0.05 GeV/c 2 



are rejected to suppress the background from J/tp 
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FIG. 1: (a) Scatter plot of M„+ n -„o versus M 7T + 7r -. (b) 
Scatter plot of M,o,o,o versus M^o w o (3 entries per event). 

For J/tp —> 73-7T , the candidate events are required 
to have no charged track. The ir° — > 77 candidates 
are formed from pairs of photon candidates that are 
kinematically fit to the tt° mass, and the \ 2 horn the 
kinematic fit with 1 degree of freedom are required to 
be less than 25. True ir° mesons decay isotropically 
in the 7r° rest frame, and their decay distributions 
are flat, contrary to 7r° candidates originating from 
wrong photon combinations. To remove wrong photon 
combinations, the decay angle, defined as the polar 
angle of a photon in the 7r° rest frame, is required to 
satisfy \cos6d ecay \ < 0.95. Events with at least seven 



and less than nine photons, which form at least three 
distinct tt° candidates, are selected. A 7C kinematic fit 
is performed to the J/tp — > 737r° hypothesis (constraints 
are the 4-momentum of J/tp and the three 7r° masses), 
and X7C < 60 is required. If there is more than one 
combination, the combination with the smallest X7C * s 
retained. Events with |M 7 „.o — to w | < 0.05 GeV/c 2 are 
rejected to suppress the background from 3/tp — > unr°ir°. 

The distributions of the selected events in the 
A/ 7r + 7r - 7r o-M 7r + w - and M^o^o^o-M^o^o planes are shown 
in Fig. [T] (a) and (b), respectively. The clusters 
corresponding to 77 — > 3ir, rj' — > 3tt and 7/(1405) — > 
/o(980)7T° can be clearly discerned; the lu signal is 
also evident in Fig. Q] (a), which mainly comes from 
the background channel J/tp — >• W7r°, while it is not 
observed in the neutral channel, as expected from charge 
conjugation symmetry. 

To confirm that the apparent signal for ry(1405) — > 
/o(980)7T° is not caused by background events, we 
perform a study with an inclusive MC sample of 2.25 x 10 s 
J/ip events generated according to the Lund-Charm 
model 23 and the Particle Data Group (PDG) decay 
tables 26j. After the same event selection as above, 
neither ?y(1405) nor /o(980) are seen in the mass spectra. 
Non-/o(980) or non-77(1405) processes are studied using 
the / (980) sidebands (0.88 GcV/c 2 < M 7r+w - (w o 7r o ) < 

-(7r°7r°) < 



0.93 GeV/c 2 and 1.05 GeV/c 2 < M„ 



71- 777 channel must be larger than that of the J/?/' M 



1.10 GcV/c 2 ) or the 7/(1405) sidebands (1.15 GeV/c 2 < 



+7, 



< 1.25 GeV/c 2 and 1.55 GeV/c 2 < 



M 7r + 7r - 7r o( 7r o, r o 7r o) < 1.65GeV/c ). No peaking structures 
are observed. 

Fig. [2] (a) and (b) show the 7r + 7r~ and 7r°7r° 
mass spectra with the requirement 1.3 GeV/c 2 < 
7l/ 7r + 7r - 7I .o( 37I .o) < 1.5 GeV/c 2 . A fit is performed 
to the 7r + 7r~ mass spectrum with the /o(980) signal 
parameterized by a Breit-Wigner function convolved with 
a Gaussian mass resolution function plus a second-order 
Chebychev polynomial background function. The mass, 
width and number of events of the /o(980) obtained from 
the fit are m = 989.9±0.4 McV/c 2 , T = 9.5±1.1 McV/c 2 
and N = 706 ± 41, respectively. A fit to the 7r 7r° mass 
spectrum, shown in Fig. [2] (b), is performed in a similar 
fashion. The mass, width and number of events of the 
/o(980) obtained from the fit are to = 987.0±1.4 MeV/c 2 , 
r = 4.6±5.1 MeV/c 2 (less than 11.8MeV/c 2 at 90% C.L.) 
and N = 190 ± 30, respectively. The measured width of 
the /o(980) is much narrower than the world average. 

Figs. [3] (a) and (b) show the tt + tt~tt q and 7r 7r°7r° 
mass spectra where 7r + 7r _ (7r°7r°) is in the /o(980) mass 
region (0.94 GeV/c 2 < Af I+I - (l o #) < 1.04 GeV/c 2 ). 
In addition to the ?y(1405), there is an enhancement 
at 1.3 GeV/c 2 . A fit is performed to the 7r + 7r~7r° 
mass spectrum. The two peaks are parameterized 
by efficiency-corrected Breit-Wigner functions convolved 
with a Gaussian resolution function. The mass and 
width of the small enhancement are fixed to PDG 
values of /i(1285) [2q |. The background is described 
by a third-order Chebychev polynomial with shape 
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FIG. 2: The tt n~ and n tt invariant mass spectra with 
7T + 7r _ 7r (37r°) in the 77(1405) mass region. The solid curve is 
the result of the fit described in the text. The dotted curve is 
the /o (980) signal. The dashed curve denotes the background 
polynomial. 

parameters determined from a simultaneous fit to the 
7r + 7T~7r° mass spectrum while the 7r + 7r _ invariant mass 
is found in the /o(980) sidebands. The normalization of 
each component is allowed to float. The masses, widths, 
number of events, efficiencies, and the product branching 
ratios of the 7/(1405) and the possible /i(1285)/t)(1295) 
contribution arc listed in Table |U The statistical 
significance of the ?7(1405) is determined by the change 
of the fit likelihood -2\nL obtained from the fits with 
and without the assumption of the ?/(1405) and is found 
to be well above 10tr. The significance of the potential 
/i (1285) / 7/(1295) contribution is determined to be 3.7er. 
A fit to the Tr°ir°ir mass spectrum is performed in a 
similar fashion, shown in Fig. [3] (b). The significance of 
the 7/(1405) is determined to be larger than IOct. For a 
possible /i(1285)/r/(1295) contribution, the significance 
is only 1.2cr, and we derive an upper limit on the 
branching ratio at the 90% C.L. using the Bayesian 
method. 

An angular-distribution analysis is performed with the 
selected J/ip -> 77/(1405) -> 7/ (980)7r° — » 77r+7r-7r° 
events. Backgrounds arc subtracted using the /o(980) 
sideband events. For radiative J/ip decays to a J p = 0~ 
meson, the polar angle 9 1 of the photon in the J/ip 
center-of-mass system should be distributed according 
to 1 + cos 2 # 7 . In the case of a J p = 1 + meson, the 
distributions d ^ sQ ~ 1 + 2|a| 2 + (1 - 2\a\ 2 ) cos 2 6> 7 and 
dcoJ; o(980) ~ 2 + (M 2 - 2)sin 2 /o(98O) are expected, 
where #/ (98o) is the polar angle of /o(980) in the helicity 
frame of 7/(1405), a is the ratio of helicity 1 to helicity 
0. For the J p = 1 + assumption, \a\ 2 is determined to 
be 2.10 ± 0.26 from a fit to cos6» /o(980) (Fig. [3] (c)). The 
fitting x 2 /n.d.f. of Fig.[3](d) with |a| 2 = 2.10 is 59.4/15. 
For the J p = - assumption, the fitting -y^/n.d.f. of 
Fig. [3] (d) is 38.4/15. Comparing the probability of 1 + 
hypothesis to 0~ hypothesis, the ratio of the probabilities 
is 4.1 x 10~ 4 . The fitting results favor the J p = 0" 
assignment of the 7/(1405). 

Figs. U] (a) and (b) show the 7r + 7r~7r° and 7r 7r°7r° mass 
spectra in the 77' region. A fit is performed to the 7r + 7r~7r° 
mass spectrum, shown in Fig. [J] (a). The shape of the 7/' 
is obtained from MC simulation and the mass and width 
of the 7/' are fixed to its PDG values j2(|. The shape 




M(f (980)ii°)(QeV/c 2 ) ' M(f (980)n°)((3eV/c 2 ) 




FIG. 3: Results of the fit to (a) the / O (980)(7r + 7r")7r° and 
(b) /o(980)(7r°7r )7r° invariant mass spectra. The solid curve 
is the result of the fit described in the text. The dotted 
curve is the /1 (1285) /r/ (1295) and 7/(1405) signal. The 
dashed curves denote the background polynomial. Angular 
distributions of the signal include efficiency corrections, (c) 
The cos Of (9 8 o) distribution. The fitting result of cos#j- (g g0 ) 
is |a| 2 = 2.10 ±0.26. (d) The cos0 7 distribution. The solid 
line is the prediction for the J p = 0~ hypothesis, and the 
dashed line is the prediction for the J p = 1 + hypothesis with 
|q| 2 = 2.10. 



of the peaking backgrounds [J/ip — > 77/' — > 77 p° (n + ir~ ) 
and J/ip — > 77/' — > 77w(7r + 7r~7r )) are from exclusive 
MC samples with predicted background levels of 361±32 
and 32±6 events, normalized by branching ratios in the 
PDG [2(| and fixed in the fit. The error on the number of 
events is estimated by changing the normalization by one 
standard deviation from the PDG value. A second-order 
Chebychev polynomial is used to describe the sum of 
other non-peaking backgrounds. For 7/ — > 7r 7r°7r°, 
X 2 (77r 7r°7r°) < X 2 (7W°tt°) and |M 77 - m v \ > 0.03 
GeV/c 2 are additionally required to remove background 
events from 7/' — > r)ir Q ir°. A fit to the 7r 7r°7r° mass 
spectrum is performed as shown in Fig.[4](b). The shapes 
of the 7/' and non-peaking backgrounds are described 
analogously to the charged mode. The efficiencies and 
the product branching ratios for the 7/' obtained from 
the fit are also listed in Table |TJ From our measurement 
and the world average values for branching ratios of 
if -> t/tttt HI, we determine r± = (8.87 ± 0.98) x 10~ 3 
and r Q = (16.41 ± 1.94) x 10~ 3 . 

The systematic uncertainties on the signal yield arise 
from fit ranges, signal shapes and background estimation. 
In detail, for the 7/(1405) signal, the signal shape is 
given by a Breit-Wigner function with floating mass 
and width. Its uncertainties are estimated by fixing 
the mass and width of the 7/(1405) to the world 
average values. For the possible /i(1285)/t/(1295) 
contribution and 7/' signal, uncertainties arc estimated 
by changing their mass and width by one standard 
deviation from the PDG values. The uncertainty due 
to the assumed background shape for the 7/(1405) and 
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TABLE I: Summary of measurements of the masses, widths, number of events, the MC efnciencies(e) and the product branching 
ratios of Br (J/ip -¥ -yX)xBr(X ->■ 7r°/o(980))xBr(/ (980) ->■ tvtt) and the decay branching ratios of Br(Y -¥ 3vr) for the 
7T + 7r~7r° and 7r 7r°7r channels, where X represents 77(1405) and the possible fi (1285) /t) (1295) contribution, Y represents t] . 
Here for the branching ratios, the first errors are statistical and the second ones are systematic. 



Resonance M(McV/c^) r^MoV/V) N event t(%) Branching ratios 

77(1405)(7r+7r-7r u ) 1409.0 ±1.7 48.3 ± 5.2 743 ± 56 22.20 ± 0.21 (1.50 ± 0.11 ± 0.11) X lO" 5 

r;(1405)(7r 7r 7r ) 1407.0 ±3.5 55.0 ±11.0 198 ± 23 12.83 ±0.11 (7.10 ± 0.82 ± 0.72) X 10" 6 

/l(1285)/77(1295)(7r+7r-7r°) fixed fixed 60 ± 18 26.99 ± 0.23 (9.99 ± 3.00 ± 1.03) X 10" 7 

/l(1285)/r?(1295)(7r 7r°7r°) fixed fixed 23 16.75 ±0.13 < 7.11 X 10" 7 

ri'(ir+TT-ir ) fixed fixed 1014 ± 39 22.52 ±0.19 (3.83 ± 0.15 ± 0.39) X 10" 3 

V(7r O 7r 7r°) fixed fixed 309 ± 19 7.57 ± 0.08 (3.56 ± 0.22 ± 0.34) X 10~ 3 
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FIG. 4: Results of the fit to (a) the TT + n~n° and (b) 3tt° 
invariant mass spectra. The solid curve is the result of the 
fit described in the text. The dotted curve is the 77' signal. 
The dashed curves denote the background polynomial. The 
dash-dotted curve in (a) describes the peaking background. 



the /i (1285) / ?7(1295) has been estimated using different 
Jo (980) sidebands; while that of the 77' is studied using 
different order polynomials. For rf — > 7r + 7r~7r°, the 
uncertainty of the peaking background is estimated by 
using the shape from the ir° sidebands instead of using 
the shape from exclusive MC samples. The systematic 
errors on the branching ratio measurements are also 
subject to systematic uncertainties in the number of 
J/ip events [2(|, the intermediate branching ratios [26| . 
the data-MC difference in the it tracking efficiency, 
the photon detection efficiency, particle identification, 
the kinematic fit and the ir° selection. Combined in 
quadrature with the uncertainty from the mass spectrum 
fitting, the systematic errors on the product branching 
ratios of the 77(1405), /i(1285)/r?(1295) and 77' are 
summarized in Table [T] 

The ?7(1405) -> /o(980)7T° signal could arise via 
77(1405) -> a%(980)n° and ag(980) - / (980) mixing. 
Using the branching ratio of J/ip —> 777(1405) — > 
7777r + 7r~ and the largest PDG value of T(77(1405) — > 
a (980)7r)/r(77(1405) -> tjtttt), Br(J/ip -> 777(1405) -)• 
7ag(980)7r° -> 7 7r ?77r ) = (8.40 ± 1.75) x 10~ 5 . The 
ag(980) - / (980) mixing intensity (£ o/ = Br( Xc i -> 
/o(980)7T° -> tt+tt-tt )/ Br (xd -> a°(980)7T° -> rm°TT )) 
measured at BESIII is less than 1% (90% C.L.) The 
branching ratio of J/ip — > 7?7(1405) — > 7ao(980)7r° — > 
7/o(980)7T° — > jir + ir~ir° is thus expected to be less than 
(8.40±1.75) x 10 -7 , which is much smaller than the result 
that we measure. Therefore, ag(980) — /o(980) mixing 
alone can not explain the observed branching ratio of 



7?(1405) -> / (980)^°. 

To interpret the anomalously narrow width of /o (980) 
and large isospin violation observed in this channel, 
J.J. Wu et al. recently propose that a novel scenario 
of triangle sing ularity would play a crucial role in 
this process [29j. Further theoretical and experimental 
studies are needed for a better understanding of the 
underlying dynamics. 

In summary, we have studied J/ip — > "f3n decays. 
The isospin violating decay 77(1405) — > /o(980)7r° is 
observed for the first time with a statistical significance 
larger than 10er in both the charged and neutral 
modes. According to our measurement, the ratio of 
5r(?7(1405) -> / (980)7r° ->■ tt+tt-tt ) to Brhj{U05) -> 
ag(980)7r° -)• ?77r 7r ) is (17.9 ± 4.2)% [H, which is 
one order of mag nitude larger than the a[J(980) - / (980) 
mixing intensity (less than 1%) determined at BESIII 
previously (28|. The measured width of the /o(980) is 
anomalously narrower than the world average. There is 
also evidence for an enhancement at around 1.3 GeV/c 2 
(potentially from the /i(1285)/t7(1295)) seen with a 
significance of 3.7ct in the charged mode and 1.2cr in 
the neutral mode. For the decay 77' — > ir + ir~ir°, the 
branching ratio that we measure is consistent with the 
CLEO-c measurement [30| , and the precision is improved 
by a factor of four. For the decay 77' — > it is two times 
larger than the world average value [26| . Using our new 
measurement of the decay rates for 77' — > 37r, the values of 
r± and ro are more than four standard deviations away 
from both -n -0 -?? mixin g pr ediction and the chiral unitary 
framework prediction [19( . 
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